Localization of correlated fermions in optical lattices with speckle disorder 
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Strongly correlated fermions in three- and two-dimensional optical lattices with experimentally 
realistic speckle disorder are investigated. We extend and apply the statistical dynamical mean- field 
theory, which treats local correlations non-perturbatively, to incorporate on-site and hopping-type 
randomness on equal footing. Localization due to disorder is detected via the probability distribution 
function of the local density of states. We obtain a complete paramagnetic ground state phase 
diagram for experimentally realistic parameters and find a strong suppression of the correlation- 
induced metal insulator transition due to disorder. Our results indicate that the Anderson-Mott 
and the Mott insulator are not continuously connected due to the specific character of speckle 
disorder. Furthermore, we discuss the effect of finite temperature on the single-particle spectral 
function. 
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I. INTRODUCTION 

Anderson localization^^ as well as strong correlation 
effects due to interactions^ have been studied intensely 
for decades in solid state physics. The elaborate inter- 
play between these two fundamental phenomena gives 
rise to challenging questions'^' for both experiment and 
theory. Recently, experiments with ultracold atoms in 
disordered lattices have opened an alternative route of 
investigation to ga in a better understanding of this fun- 
damental interplay.™^ A high degree of tunability of the 
disorder and interaction strength is inherent to these ex- 
periments. 

Three different ways of realizing disordered optical 
lattices have been proposed: i) by loading an addi- 
tional atom ic sp ecies into the lattice and suppressing 
its mobility! 11 * 12 * Thereby a binary disorder potential 
is realized due to the interatomic interaction; ii) by 
superimposing two laser beams with incommensurate 
frequencies j 1 3 1 14 l thereby generating a stationary quasi- 
periodic lattice; iii) by using an op tical speckle fielcH^ 
superimposed onto an optical lattice P^EZI 

In general, particles in a disordered lattice are ex- 
posed to an on-site potential varying from site to site, 
which already gives rise to localization due to coherent 
backscattering. 1 In tight-binding lattice models this ran- 
domness is termed diagonal disorder. Additionally in 
a tight-binding description, the hopping amplitude be- 
tween two lattice sites depends on the specific disorder 
potential as well, giving rise to off-diagonal disorder 
A complete description of particles in a disordered lattice 
should incorporate both diagonal and off-diagonal disor- 
der, i.e. random on-site energies and hopping disorder 
already in a single particle picture. 

A realistic theoretical description of ultracold atoms 
in a speckle disordered optical lattice should be capable 
of describing interaction and disorder on equal footing, 



preferably non-perturbatively. Furthermore, it should in- 
corporate diagonal and off-diagonal disorder and, last but 
not least, should be able to treat other experimental fea- 
tures, such as finite temperature and the presence of the 
harmonic trap. In the c ase o f bosons, there have been 
several theoretical works 24 25 combining some of these 
requirements. To the best of our knowledge, comparable 
studies have so far not been performed for fermions. 

The aim of our paper is to close this gap and to 
provide a realistic description of strongly correlated 
fermions in a disordered optical lattice. For this purpose, 
we empl oy t he statistical dynamical mean-field theory 
(DMFT)P^H28] tQ golve the Anderson-Hubbard Hamil- 
tonian numerically. The statistical DMFT incorporates 
both strong correlations and disorder-induced fluctua- 
tions and is known to give accurate results for high- 
dimensional systems, i.e. d > 3. We extend the statisti- 
cal DMFT scheme to take hopping disorder into account. 
Realistic distributions for the Hubbard parameters, cal- 
culated by S. Q. Zhou and D. M. Ceperley™ are used to 
describe the speckle disorder potential. 

As our main result we find that the paramagnetic 
ground state phase diagram differs strongly from the 
phase diagram for pure diagonal box (homogeneous) dis- 
order, as obtained within typical medium theory ! 30 " 33 * In 
particular, the Mott insulator and the Anderson-Mott in- 
sulator are not continuously connected in case of speckle 
disorder. Our results are reappraised by a complemen- 
tary investigation of the two-dimensional square lattice 
by means of real-space DMFT (RDMFT) ISSlMlMJ 



II. MODEL OF CORRELATED FERMIONS IN 
SPECKLE DISORDERED OPTICAL LATTICES 

Ultracold fermionic atoms, such as 6 Li or 40 K, in dis- 
ordered lattices are described by the Anderson-Hubbard 
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Hamiltonian 



-E 



UiTi^n^ (1) 



where c\ a (q ct ) denotes the creation (annihilation) oper- 
ator at a lattice site i with spin a = ±1/2. The fermionic 
number operator is given by = c\ a Ci a . The hopping 
amplitude between sites i and j is denoted by tij, the 
local interaction potential is parametrized by Ui, and the 
chemical potential is given by \i. In the following, we con- 
sider fermions on a bipartite lattice with a semi-elliptical 
model density of states (DOS) pSHEH which is character- 
ized by the connectivity K. It is related to the lattice 
coordination number Z via K = Z — 1. The hopping 
amplitude Uj is assumed to be only non-zero between 
nearest neighbor sites i and j. 

In experiments with cold atoms, the speckle disorder 
potential is created by a co herent la ser beam that is scat- 
tered by a diffusor plate f^E 17|39 | 4Q| A statistical analysis 
of the scattering process^ shows that the probability dis- 
tribution function (PDF) of the resulting light intensity 
pattern obeys p 7 (J) = 6(7) exp(-J/ (/))/(/), where (I) 
is the averaged light intensity and Q(x) denotes the Heav- 
iside function. By superimposing the speckle light pat- 
tern onto the optical lattice, the atoms are subjected to 
a random optical dipole potential Vo(r) ex I(r)p^ which 
is attractive for red-detuned laser light or repulsive for 
blue-detuned laser light. We will consider the latter, i.e. 
a repulsive potential. Within the tight-binding model 
description this random potential gives rise to diagonal 
disorder, i.e. random on-site energies e^, which are drawn 
from a PDF p e {ti) given by 
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ex P(-x)°( e )> 



(2) 



where A denotes the disorder strength. We note that this 
PDF of the on-site energies is unbounded from above, 
in contrast to other typically used distributions, such as 
box or binary disorder. Furthermore, we assume that the 
on-site energies of all lattice sites are independently and 
identically distributed. Within a tight binding model, 
the speckle disorder potential leads to off-diagonal disor- 
der as well! 16|29 lThis so-called hopping disorder manifests 
itself in the random hopping amplitudes . For a given 
disorder distribution, inducing fluctuations in both the 
hopping and the on-site energies, the hopping coefficient 
Uj at a neighboring pair of sites is correlated with the 
difference in on-site energies Ae = — ej and a realistic 
description requires modeling using a joint PDF 2 ^ 



p Ae ,t(Ae,£) ^PAe(Ae) ■ p t (t) . 



(3) 



In our model ([!]), the joint PDF £>Ae,t(Ae, t) is incor- 
porated based on the data by S. Q. Zhou and D. M. 
CeperleyP^ We account for the dependence of these two 
random variables by a conditional PDF p t (t\/S.e) for the 
hopping t after the on-site energies on the respective sites 
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FIG. 1: (Color online) One-dimensional illustration of the 
disordered lattice problem. The speckle field induces random 
on-site energies given by the probability distribution func- 
tion p e (e). The hopping amplitudes Uj are random variables, 
which manifests itself in the probability distribution function 
PAe,t(Ae, £). The local interaction strength U% depends on the 
local Wannier function, which results in a further probability 
distribution function p e ,u(£, U). 



were sampled. The details are given within the introduc- 
tion of the method in subsection |III C| and the Appendix 
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Furthermore, for a short-range interaction between 
particles the local Hubbard interaction parameter is pro- 
portional to the integral over the fourth power of the 
Wannier function on the same site, which in turn depends 
on the random lattice potential. Hence, the on-site in- 
teraction coefficient U is a random variable as well, and 
the joint PDF p ej u(e,U) of the on-site interaction and 
the on-site energies also needs to be taken into account. 

Similar to the joint PDF of the difference in the on- 
site energies and the hopping coefficients, the joint PDF 
p e? c/(e, U) is incorporated based on data from Ref. l29land 
accounted for by the conditional PDF pu(U\e). The me- 
thodical details are explained in Sec. |III C| and Appendix 
lAl A one-dimensional illustration of the disordered lat- 
tice problem considered here is shown in Fig. [I] 

In experimental realizations, the disorder strength A is 
proportional to the speckle field strength sd = (Vo(r))P3 
and can therefore be tuned. The proportionality con- 
stant scales monotonically with the ratio of the speckle 
field autocorrelation length to the typical spatial extend 
of the Wannier function, and therefore depends on the 
experimental optical setup. In the experiment by White 
et a/., 16 the relation A = 0.97s D was found. 29 The mean 
value of the interaction strength U can be tuned by ad- 
justing the s- wave scattering length a between the two 
fermionic components by a Feshbach resonanc e 1 42 * 43 * or 
by varying the dimensionless lattice depth sl, which in 
turn also influences the hopping amplitude. Therefore, 
keeping the lattice depth fixed, the disorder strength 
and the interaction strength can be tuned independently 
by varying the speckle field strength and the magnetic 
field. The relevant energy scale for the lattice depth sl 
and the speckle field strength S£>, is the recoil energy 
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Er = hk 2 /2m, where k is the wave vector of the optical 
lattice and m is the atomic mass. Throughout this paper 
we work in energy units of the noninteracting bandwidth 
Wo = ^t\fK of the homogeneous system, where the hop- 
ping amplitude t can be related to the lattice depth si by 
an exact band structure calculation, as for instance per- 
formed by M. Greiner!^ The s- wave scattering length a 
is given in units of the Bohr radius clq. 



III. METHODS FOR SOLVING THE MODEL 

Two conceptual approaches using DMFT to disordered 
systems, both referred to as statistical DMFT, were orig- 
inally introduced in Ref. [26j In our opinion a distinction 
in terminology is justified and useful: Firstly, RDMFT, 
leading to a set of self-consistency equations for a fixed 
disorder realization or any inhomogeneity, thus consti- 
tuting a deterministic approach. This method is appli- 
cable to any finite lattice structure and incorporates dis- 
order or inhomogeneities non-perturbatively. Secondly, 
we refer to statistical DMFT as the intrinsically statisti- 
cal approach to disordered systems on the infinite Bethe 
lattice. Here, the disordered system is investigated on a 
fully stochastic level and the self-consistency applies on 
a level of PDFs for the Green's function. The statistical 
DMFT as well as the RDMFT are described in Sees. UTTAl 



and IIIB In Sec. Ill C we explain how we incorporate 
the joint PDFs of the Anderson-Hubbard model. 



A. Statistical dynamical mean-field theory 

The statistical DMFT^ is a self-consistent computa- 
tional scheme for determining the PDF of the local single- 
particle Green's functions, i.e. p [Ga a (jjj)]. Here, Gij a (uo) 
is the Fourier transform of the retarded Green's function 
Gija(t) = -i6(t)([c icr (t),ct o .(0)] + ), where [.., ..] + denotes 
anticommutator brackets. In the following, only para- 
magnetic solutions of the Anderson-Hubbard model are 
considered. Therefore the spin index a is omitted here- 
after. 

In the absence of interactions, the renormalized per- 
turbation theor}El shows that the local Green's function 
can always be expressed as 



Ga(ui) 
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Ti{uj)+iri 



(4) 



where the hybridization function Ti(uS) describes the cou- 
pling of site i to nearest neighbor lattice sites. For nu- 
merical reasons we introduce a finite broadening factor 
r] > 0. The limit 77 — > is important to detect local- 
ization. An additional self-energy ^^(u) accounts for 
interaction effects and is defined via a Dyson equation. 
In the dynamical mean-field theory^ approximation the 
self-energy is local, i.e. S^(o;) — » E(o;)^j, and the local 



single-particle Green's function given by 

1 

Gu(u) = — ■ 



(5) 



The local approximation for the self-energy becomes ex- 
act in infinite dimensions, as was shown by W. Metzner 
and D. VollhardtPSl and is known to be a good approxima- 
tion in three spatial dimensions. In a disordered system, 
translational invariance is broken and the local single- 
particle Green's function varies randomly from site to 
site, giving rise to a PDF p[Ga(uj)]. Within statistical 
DMFT, this PDF is determined by an ensemble with a 
large number N of Green's functions. 

On the Bethe lattice, the hybridization function is 
given as sum over diagonal cavity Green's functions 
gCO^ BBBZBSHZ I^.e. 



(6) 



The corresponding cavity hybridization functions (u) 
are given as a sum over z— 1 = K cavity Green's functions 
G§ j \u). In this case, the exclusion of site i is no longer 
relevant, since the Bethe lattice is loop-free. In a similar 
manner, subsequent equations reproduce the structural 
dependence of a sum over K diagonal cavity Green's 
functions. An approximation first successfully applied 
by R. Abou-Chacra et aZ.p^and used in other works^M^ 
is that the cavity hybridization function is given by the 
hybridization function of the initial lattice 
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with the sum now extending over K Green's functions. 

The local self-energy is determined by the solu- 

tion of the corresponding impurity problem, defined by 
the given hybridization I\ (lj) pOSO Therefore, the statis- 
tical DMFT maps the original lattice model onto an en- 
semble of impurities, whose PDF has to be determined 
self-consistently. 

In practice, the self-consistent calculation scheme con- 
sists of the following steps: The starting point is an ini- 
tial PDF p[Gu(jjj)] and the calculation is performed by 
(i) drawing a random on-site energy from the PDF 
p e {€i) given in Eq. ([2| and a random on-site interaction 
Ui from the conditional PDF pu(U\e) for each sample of 
the ensemble; (ii) for each sample i random hopping am- 
plitudes tij are drawn from the conditional PDF p t (t\Ae) 
depending on the difference Ae = ei — ej of the on- 
site energies and ej of the randomly determined near- 
est neighbor j; (hi) The hybridization function r^(cj), 
with the local single-particle Green's function Gjj(u) of 
the nearest neighbors randomly sampled from the PDF 
p [Ga(cj)], is determined via Eq. ^ for each sample; (iv) 
the local self-energy E^(o;) is calculated from the solution 
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FIG. 2: (Color online) Illustration of the statistical dynami- 
cal mean-field theory employed in this work. The many-body 
problem with disorder is mapped onto an ensemble of sin- 
gle impurities, which are coupled to a random ensemble of 
bath Green's functions, which is determined self-consistently. 
G n represents the nth sample from the ensemble of Green's 
functions. 



of the local impurity problem by using an impurity solver; 
(v) the local single-particle Green's function Gu(uj) is 
calculated using Eq. (pi); (vi) having calculated a com- 
pletely new ensemble {Gu(uj)} : a new PDF p[Gu(uj)} is 
obtained by construction of a histogram and we return to 
step (i). The algorithm is repeated until self-consistency 
for p[Gu(u)] is achieved. We note that this method in- 
corporates spatial fluctuations, i.e. quantum interference 
effects, caused by disorder via step (iii). Schematically 
the computational procedure is presented in Fig. [2] 

The relevant physical observable is the local density of 
states (LDOS), given by pi(uo) = — ^Im(G^(cj)), which 
is a random quantity in disordered systems. The corre- 
sponding distribution p[pi(uj)] is obtained by counting all 
values of the LDOS for each frequency and by construct- 
ing a histogram. Statistical fluctuations are minimized by 
an artificial increase of the size of the ensemble after hav- 
ing reached self-consistency. This is done by constructing 
the histogram on the basis of typically 50 — 100 succes- 
sive DMFT iterations. From this PDF the expectation 
value can be determined, i.e. the arithmetically averaged 
LDOS, 



<p(")Uh = / dp\uj)p[p\uj)}p\u 



(8) 



and the typical value, which as in typical medi um theo ry, 
is approximated by the geometrical averagd ^ 33 * 50 * 51 ! 

(pM)geo m = exp J dp'(u) p[p'(uj)} lnp'(u), (9) 



where the dependence on uj is to be understood paramet- 
rically. In the following, the cumulative PDFs 



P\p(w)] 



p\p'{u>)]dp'(u>) 



(10) 



will also be useful to characterize the disordered system. 

We close the description of the method by a short dis- 
cussion of our impurity solver. The most time-consuming 
part in the scheme is the solution of a large number of im- 
purity problems, requiring a fast impurity solver. Here, 
we use iterative perturbation theory (IPT)p"Ml w fiich 
properly reproduces the non-interacting and atomic lim- 
its and was shown to qualitatively describe the Mott- 
Hubbard metal-insulator transition at a finite critical 
interaction strength U C F** Within IPT, the self-energy 
is calculated to the second order in U in the non- 
renormalized perturbation expansion. Using this impu- 
rity solver, Green's function ensemble sizes of the order 
N ~ 10 3 are computationally feasible within a paral- 
lelized code. 

The original IPT, which was restricted to the half- filled 
case, was subsequently extended to densities away from 
half-filling. It is now comm only referred to as modified 
perturbation theory (MPT)WMl The self-energy within 
MPT is given bj^L 



E(cj) = Un 
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(11) 



where ^^(uj) is the second order perturbation con- 
tribution to the self-energy!^ Within the perturbation 
expansion, the non-renormalized Hartree-Fock Green's 
function^ 



U{ni) -Ti(Lo)+ir} 



(12) 



are used as propagators. The parameter p is fixed by 
requiring (m) = (ni) nF W^ The parameters a and b are 
determined such that the first three spectral moments 



uj m p(u)duj 



(13) 



with m = 0, 1, 2 are reproduced and that the atomic limit 
is recovered correctly.^ For further details, the reader is 
referred to the work by M. Potthoff et al 52 and our earlier 
work. 28 



B. Real-space dynamical mean-field theory 

Fermions with a semi-elliptic DOS in many ways ex- 
hibit qualitatively the same physics as fermions in three- 
dimensional lattices. It is however important to note 
that qualitative differences may arise in lower dimensions. 
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FIG. 3: (Color online) Color coded probability distribution 
function (histogram) p e ,u(c, U) of the on-site energies and the 
on-site interaction for speckle field strength sd = 0.05Er and 
s-wave scattering length a = lOOao for a 3d lattice. The 
lattice depth is given by sl = IOEr. 



For this reason we also investigate fermions on a square 
lattice within RDMFT.^ 27 ' 34 ^ 5 ' Besides describing the 
Mott-Hubbard metal-insulator transition and magnetic 
order, RDMFT is also capable of treating spatial inho- 
mogeneities such as disorder. As in DMFT, each lattice 
site is mapped onto a single-impurity Anderson Hamil- 
tonian within RDMFT, where the hybridization function 
has to be determined self-consistently. 

Starting with an arbitrary hybridization function, the 
solution of each impurity problem is provided by MPT, 
as in statistical DMFT, and leads to a set of on-site self- 
energy functions T, l Sij(iuj n ). They determine the self- 
energy matrix in the real-space representation 

(S)« = . (14) 



C. Joint probability distribution functions for the 
Hubbard parameters 

The Hubbard parameters used both in the statistical 
DMFT and the real-space DMFT simulations are ran- 
domly generated based on the distributions obtained by 
S. Q. Zhou and D. M. Ceperley^ within an imaginary 
time evolution approach for the Wannier functions. Al- 
though these calculations were originally performed for 
bosonic 87 Rb, the random parameters e^, j depend 
only on the structure of the single particle states and can 
thus also be used for fermionic system after an appropri- 
ate rescaling. 

The three sets of parameters for the Hubbard model 
underly statistical fluctuations. For the 3D case with 
equal laser intensity along each of the three axes, all pa- 
rameters are unique functions of sl and sd for a given 
atomic species, when expressed in units of Er. This is 
not the case in the two-dimensional lattice, relevant for 
our real-space DMFT calculations. Here, the interaction 
parameter U depends on the shape and strength of the 
axial trapping potential, which may vary significantly in 
different experiments. In this anisotropic case with fixed 
a s , the lattice depth sl does not uniquely characterize 
the point in the phase diagram. Here we therefore give 
energies in units of the noninteracting bandwidth Wq. 

The on-site energies and interaction parameters Ui 
are sampled from the data given in Ref. 29 inherently 
containing correlations, as shown in Fig. [3] Up to a very 
good approximation, the on-site energies and tunnel- 
ing parameters tij are independent, however there is a 
significant correlation between the difference in on-site 
energies Ae^ = — ej and the tunneling parameter , 
as is shown in Fig. [4] These correlations are taken into ac- 
count both within the statistical DMFT and the RDMFT 
simulations by sampling t from a conditional distribution 
p(t|Ae) for a given (previously sampled) Ae. More tech- 
nical details of how the conditional PDFs are constructed 
are given in Appendix [Aj 



Following the Dyson equation, the interacting lattice 
Green's function is given by 



G(iLj n ) 1 = G (iuj n ) 1 - T,(iu n ) 



(15) 



where uo n are the Matsubara frequencies. The non- 
interacting Green's function Go(iuj n ) in real-space rep- 
resentation is given by 



G (iu n ) 



(/i + icj n )l - t - V, 



(16) 



where 1 is the unity matrix, t is the matrix of hopping 
amplitudes, and (V)ij = £i$ij denotes the matrix of ran- 
dom on-site energies. Together with E q. ([5] ) and the di- 
agonal elements from inverted relation (15|~a set of local 
hybridization functions T^(uj) is extracted, which closes 
the self-consistency loop. 



IV. RESULTS 

In our calculations we consider ultracold 40 K atoms 
in a mixture of two hyperfme states in an optical lattice 
generated by a laser with a wavelength Xl = 738nm. 
The lattice depth is fixed to = WEr. Moreover, 
we consider the half-filled case, i.e. band filling v = 
Jr Xli C r( n * cr ) = 1> which is accomplished by adjusting the 
chemical potential \i. For the statistical DMFT calcula- 
tions, the lattice connectivity K = 6 was chosen. 

Describing the physics of strongly correlated fermions 
in a speckle-disordered optical lattice necessitates a 
proper definition of the relevant phases. The Mott in- 
sulator is incompressible and its spectrum is gapped. Its 
correlation gap is proportional to the interaction strength 
U. A gapped spectrum in turn means that the arith- 
metic average of the LDOS at the Fermi level vanishes 
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FIG. 4: (Color online) Color coded joint probability distri- 
bution function (histogram) pAe,t(Ae,£) of the difference in 
on-site energies between neighboring sites of a 3d lattice and 
the respective hopping amplitudes for a speckle field strength 
sd = OAEr and lattice depth sl = IOEr. As restricted 
by the symmetry of Uj under the exchange of 6i <H> €j, the 
distribution only depends on |Ae|. 



(p(u = 0)) arit h = 0. A second insulating phase which we 
find here is the Anderson insulator 1 for zero interaction. 
This phase is compressible and its spectrum is point-like, 
which can be attributed to an absence of diffusion.^ In 
the interacting case, we define a state to be localized or 
extended, if the spectrum of the single-particle Green's 
function is point-like or continuous respectively. The 
Anderson-Mott insulator is defined by localized single- 
particle excitations at the Fermi level. Recently, it was 
shown that localized single particle excitations imply a 
vanishing conductivity in case of weak interactions.^ 
The detection of Anderson localization within statisti- 
cal DMFT will be explained in the following subsection 
IV A Finally, the paramagnetic metal is compressible 
and therefore has a non- vanishing LDOS (p(u = 0)) arith 
at the Fermi level, where the states are extended. 




FIG. 5: (Color online) Comparison of arithmetically (black 
solid lines) and geometrically (red dashed lines) averaged 
spectral functions with increasing speckle field strength 
s D (s D = 0.0Er,0.1Er,0.2Er,0.3Er,0AEr) in the non- 
interacting limit a — 0. Parameters are v — 1.0, sl — 10 Er, 
and n = 10~ 3 . 



the arithmetic mean, is that the geometric mean of the 
LDOS converges to zero with increasing A. This can be 
attributed to the transition from extended states to local- 
ized stat es, which i s also observed within typical medium 
theory ffi* 31 * 33 * 50 * 5 -^ However, not only averages, but also 
the full PDF p[p(uo)] is accessible within the statistical 
DMFT, which enables a more accurate detection of lo- 
calization than attributing it to a vanishing geometrically 
averaged LDOS. This approach will be explained in de- 
tail in the next paragraph. The PDFs p[p(uj)] associated 
with the data shown in Fig. [5] are plotted in Fig. [6] As 
expected in the non-disordered case (cf. panel (a)), the 




A. Detecting Anderson localization within 
statistical DMFT 



The arithmetic and geometric averages of the LDOS 
obtained within statistical DMFT, which for the 
non-inter acting case reduces to the local distribution 
approach j 47 * 49 * exhibit different behavior with increasing 
disorder strength A. Fig. [5] displays the evolution of the 
arithmetically and geometrically averaged LDOS in the 
non-interacting case, when the speckle field strength is 
increased from OEr to OAEr. First of all, we note that 
with an increase in the speckle field strength, the spectra 
are broadened and long tails emerge at positive energies. 
A noticeable difference between the geometric mean and 





FIG. 6: (Color online) Color coded natural logarithm of PDF 
p[p] for increasing speckle field strength sd- (a) sd = O.Oi^R, 
(b) s D = OAEr, (c) s d = 0.2E R , and (d) s D = 0.3E R . 
Parameters are a s — 0, n — 10 -3 , v — 1.0, and sl — IOEr. 
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FIG. 8: (Color online) Behavior of the cumulative PDFs 
P[p(u) = 0)] in the noninteracting case as the broadening 77 
is decreased from 10 -3 to 10 -6 for (a) speckle field strength 
sd — OAEr and (b) speckle field strength sd = 0.7 Er. Pa- 
rameters are a = 0, v — 1.0, and sl = IOEr. 
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FIG. 7: (Color online) Color coded natural logarithm of PDFs 
p[p(oj)] for decreasing broadening 77 and speckle field strength 
so = OAEr (panels (a)-(d)) and sd = 0.7 Er (panels (e)- 
(h)) of non- interacting fermions. Parameters are v — 1.0 and 
s L = 10E R . 



PDF for each frequency is given by a delta function. For 
finite disorder strength, the PDFs extend over several 
orders of magnitude. 

Within statistical DMFT, extended and localized 
states are characterized by different behavior of the PDF 
p[p(uu)] in the limit of vanishing broadening 77 — »• O ! 28 * 49 * 
This procedure is motivated by the fact that the ex- 
tended states are given by a branch cut on the real axis 
of the single-particle Green's function, whereas the lo- 
calized states are characterized by a dense distribution 
of poles in the thermodynamic limit !^ Let us consider a 
single particle initially located at site on a finite lattice 
and let \ip n ) denote the complete set of single particle 
energy eigenstates on the lattice with eigenenergies E n . 
The local Greens function then takes on the forrrP^ 



G00M 



E 



fn 



Ern 



(17) 



Here, f n — (0|V ; n)(V ; n|0) denotes the overlap of the eigen- 
state n with the Wannier function on site 0. The return 
probability is given by po^o(t —> 00) = / 2 .^In case 



of an extended eigenstate, the residue f n is proportional 
to the inverse number of occupied lattice sites TV -1 and 
the return probability approaches zero in the limit of an 
infinite system. In contrast, for a localized state the spec- 
trum is given by a dense distribution of poles in the al- 
lowed energy interval on the real cj-axis. The residues 
approach a finite value, but some will dominate the sum 
in equation (17). If their values were sorted by value, 



their contribution would decrease exponentially for spa- 
tially localized states. In particular, if we introduce a 
small coupling 77 to a dissipative bath and consider a con- 
tour which encloses a small energy interval 5E, then the 
most probable value of the sum of the residues of poles 
enclosed by the contour will decrease exponentially pro- 
portional to the ratio A/5E. Accordingly, the most prob- 
able value of the imaginary part of the Green's function is 
proportional to 77.^ Hence, the LDOS will be highly frag- 
mented in case of localized eigenstates, characterized by 
dominating well-separated resonances. In contrast, the 
arithmetically averaged spectral function, i.e. the DOS 
of the system, does not exhibit this high fragmentation, 
since the spectral weight must be located somewhere in 
the lattice for every energy. Consequently, the maximum 
of the PDF of the LDOS of an extended state saturates 
at a finite value for 77 — >• 0, while an increasing amount 
of the PDF's weight shifts to zero in the case of localized 
states. 

In Fig. [7] the behavior of the PDFs at two speckle field 
strengths sp = OAEr and sp = 0.7 Er is plotted for a 
sequence of decreasing broadening 77 — ^ 0. The different 
behavior of the extended and localized states allows us 
to distinguish between them and the mobility edges can 
be identified as shown in panel (d) for sp = 0.4Er. The 
states between the mobility edges are extended, whereas 
the ones outside are localized. We observe that the states 
at the Fermi level are extended for speckle field strength 
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FIG. 9: (Color online) Paramagnetic ground state phase 
diagram of fermions in a speckle disordered optical lattice. 
Parameters are v — 1.0 and sl = IOEr. In the Munich 
experiment!^ (so far without disorder) scattering lengths up 
to a = 300ao can be achieved. 



sjj = 0.4 which is therefore metallic, whereas for speckle 
field strength sjj = 0.7 Er the states at the Fermi level 
are localized. This can also be seen from the 2d-plots 
of the corresponding cumulative PDFs at the Fermi level 
for the two speckle field strengths in Fig. [8j 



B. Paramagnetic ground state phase diagram of 
the interacting system within statistical DMFT 

Now we will discuss the main result of our investiga- 
tion, i.e. the paramagnetic ground state phase diagram 
of interacting fermions exposed to speckle disorder (as a 
function of s-wave scattering length a and speckle field 
strength S£>), as displayed in Fig. [9j 

In the absence of disorder (sp = 0), a Mott metal 
insulator transition is found at intermediate interaction 
strength. For the system considered here, we found the 
critical s-wave scattering length for the Mott transition 
a c = 117. 5ao- In the absence of interactions (a = 0), the 
Anderson transition occurs at sd = 0.65. In a system 
with both speckle disorder and interactions, three sepa- 
rate phases exist: Mott insulator, disordered correlated 
metal, and Anderson-Mott insulator. 



In Fig. 10 the arithmetically and geometrically aver- 
aged spectral functions are given for two different speckle 
field strengths sjj: in panel (a) for sjj = 0.05Er and in 
panel (b) for sd = O.li^R. In both cases the interaction 
strength is increased from the bottom to the top. Re- 
markably, the spectral functions evolve very differently 
for the two disorder strengths. For speckle field strength 
sd = 0.05I^r a correlation-induced metal-insulator tran- 
sition takes place at finite a. For S]j = O.li^R in the 
investigated regime up to a = 350ao, no metal- insulator 
transition was found. Instead, the Kondo peak, i.e. the 
coherent low-energy excitations, are stabilized for higher 



interaction values, whereas the lower Hubbard band is 
shifted away from the Fermi level. 

This behavior is caused by the redistribution of states 
into the Mott-Hubbard gap due to disorder. For suffi- 
ciently strong disorder, the gap is closed. The unbounded 
nature of the speckle disorder at any finite A gives rise 
to states with very high energies, although their number 
is exponentially suppressed in A -1 . This means that the 
Mott transition at finite disorder strength, which is de- 
scribed here, might even be an artifact of the finite size 
N in the stochastic Green's function ensemble. If this is 
the case, it would be an intrinsic feature of any finite size 
optical lattice as well. 

To gain further insight, we plotted the arithmetically 
averaged LDOS at the Fermi level as a function of the 



s- wave scattering length a in Fig. 11 for three differ- 
ent speckle field strengths. Due to the finite disorder 
strength, the arithmetically averaged LDOS at zero in- 
teraction is reduced, which means that the Luttinger the- 
orem is not fulfilled in presence of disorder. By increasing 
the interaction strength, the metallicity is improved for 
all three disorder strengths. At low disorder strength 
(sd = 0.05^^), the Luttinger theorem is asymptotically 
fulfilled for a ~ lOOao- This is in agreement with results 
for the Anderson-Hubbard model with box disorder p3E2] 
The metallicity is suddenly reduced for stronger interac- 
tions and finally a Mott-Hubbard transition takes place 
at a c = 145ao for sp = 0.05Er. For higher disorder 
strengths (sd = O.IEr and sd = 0.2^), no Mott tran- 
sition is found for scattering lengths up to a = 350ao- 

It is interesting to compare the qualitative structure 
of the phase diagram in Fig. [9] with the counterpart of 
ultracold bosonic atoms in speckle disordered lattices. In 
the latter case, an arbitrarily weak speckle field leads to 
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FIG. 10: (Color online) Arithmetically (black solid) and ge- 
ometrically (red dashed) averaged spectral function with in- 
creasing s-wave scattering length a for speckle field strength 
sd — 0.05Er in panel (a) and for speckle field strength 
s D — O.IEr in panel (b). Parameters are v — 1.0 and 
s L = 10E R . 
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FIG. 11: (Color online) Arithmetically averaged local den- 
sity of states at the Fermi level (p(0)) ar ith for three different 
values of the speckle field strength sd- sd = 0.05Er (red 
solid line), sd = 0.2Er (green dashed line) and sd — 0.2Er 
(blue dotted line). For comparison the local density of states 
at the Fermi level of the homogeneous, non-interacting sys- 
tem po(0) = 1.265 is included. Parameters are v — 1.0 and 
SL = 10E R . 
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FIG. 12: (Color online) Arithmetically averaged spectral 
function for increasing interaction U at fixed speckle field 
strength sd — 0.05Er obtained via statistical DMFT with 
K = 4 (black solid line) and real-space DMFT (red dashed 
line). The spin-summed filling is given by v — 1.0 and the 
lattice size within the RDMFT calculations was 24 x 24. 



a vanishing of the excitation gap and the Mott insulator 
only exists in the homogeneous system without disorder, 
in contrast to the fermionic case where a Mott insula- 
tor may exist at A > 0. Furthermore, the results pre- 
sented here differ from the results obtained within typ- 
ical medium theor}ffl^ for fermions with bounded box 
disorder. Although a derealization tendency was found 
for box disorder, a correlation-induced metal-insulator 
transition takes place at intermediate disorder strengths. 
The critical interaction strength is shifted to higher val- 
ues proportional to the disorder strength. Thus, the 
Anderson-Mott insulator and the Mott insulator were 
found to be continuously connected.^ All important dif- 
ferences between the paramagnetic ground state phase 
diagrams for box disorder and speckle disorder can be 
attributed to the unbounded nature of the speckle distri- 
bution. In particular, the result that the Mott insulator 
and the Anderson-Mott insulator are not continuously 
connected is not attributed to the methodical differences 
between statistical DMFT and typical medium theory. 



C. Results within RDMFT 

In order to investigate the speckle disordered square 
lattice and to assess our results obtained within statis- 
tical DMFT for fermions with a semi-elliptical DOS, we 
performed complementary RDMFT calculations. An ex- 
emplary comparison of the arithmetically averaged spec- 
tral functions obtained by the two methods for identical 
parameters is given in Fig. [l2| In these statistical DMFT 
calculations the connectivity K = 4 is chosen to obtain 
the same bandwidth as on the square lattice. Through- 
out this section U denotes the most probable value of the 



respective marginal PDF pu(U) = / de f p e ^jj(e 1 U). 

Both methods lead to qualitatively identical results. 
The differences in the spectral functions can be traced 
back to the differences of the simulated models. The 
statistical DMFT is employed for particles with a semi- 
elliptical DOS, whereas RDMFT was applied for parti- 
cles on the square lattice. Since the kinetic energy is 
connected to the lattice structure, the observed differ- 
ences are pronounced when the kinetic energy dominates 
over the interaction energy. Consequently, deviations in 
the distributions of the spectral weights are larger for low 
and intermediate interaction strengths (Fig. 12 U = 1.0 
and U = 2.1). On the other hand, the agreement is good 
for the strongly interacting case (Fig. [l2] U = 3.3 and 
17 = 4,2). 

Interestingly, a pseudo-gap at the Fermi level in the 
LDOS on the square lattice is found within RDMFT. 
This pseudo-gap arises for intermediate and strong inter- 
actions in the presence of disorder and is stable under 
variation of the system size. The pseudo-gap anoma- 
lies, also called zero bias anomalies, are a common fea- 
ture in two-dimensional strongly correlated systems with 
disorder EilMI a pseudo-gap anomaly was for instance 
found within a quantum Monte Carlo investigation of the 
Anderson-Hubbard Hamiltonian with box disorder. 63 64 

Arithmetically and geometrically averaged spectral 
functions calculated by RDMFT for two different dis- 
order strengths, namely sd = O.OIEr and sd = 0.05^, 
are displayed in Fig. [l3j Qualitatively, the spectral func- 
tions show similar behavior as obtained within statistical 
DMFT, cf. Fig. [10] For weak disorder (s D = 0.01^), 
metallic solutions are obtained for weak interactions. 
Raising the interaction, a Mott insulating phase is found, 
analogous to the case of a homogeneous system. On the 
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contrary, for larger speckle disorder (sp = 0.05Er) the 
LDOS remains finite at the Fermi level uj = 0, even at 
strong interaction U = 4.2. 

We note that within RDMFT for U = 3.1 and U = 4.2, 
the lower Hubbard-band shows a peaked structure at low 
disorder. This feature cannot be exclusively identified 
with physical properties of the system because of numer- 
ical uncertainties. Within RDMFT, an artificial broad- 
ening r] is applied, which is scaled proportionally to the 
system size, i.e. r] ex 1/L 2 . Since the peaks are not fully 
recovered for other lattice sizes, we conclude that they 
are finite size effects. 

The real space resolution of the LDOS gives us insight 
into localization effects of the system. In Fig. 14 the 
LDOS Pi{oS) is plotted for a 24 x 24 lattice and different 
interaction strengths, each for a different disorder real- 
ization. At U = 1.0 and weak disorder sd = O.OIEr 
(Fig. 14 panel (a)) the spectral weight around the Fermi 
level uj = remains finite at each lattice site. The vast 
majority of single particle states are extended and the 



(a) s n =0.01E R , u=i.o 
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system is in the metallic phase. At U = 3.1 (Fig. 
panel (b)) Hubbard bands are formed and the spectrum 
exhibits a gap proportional to the interaction strength for 
all lattice sites, indicating that the system is in a Mott 
insulating state. However, as the speckle field strength is 
increased to sjj = OAEr, the states with spectral weight 
at the same frequency are spatially separated. In other 
words, the spectrum is highly fragmented and each lo- 
cal spectrum consists of isolated delta peaks, consistent 
with Anderson localized states in the infinitely large sys- 
tem (cf. Fig. 14 panel (c)). 
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FIG. 13: (Color online) Arithmetically (black solid) and ge- 
ometrically (red dashed) averaged spectral function for in- 
creasing interaction strength U at speckle field strength sd = 
O.OIEr in panel (a) and sd = 0.05Er in panel (b) obtained 
by RDMFT on a 24 x 24 square lattice. The spin-summed 
filling is given by v — 1.0. 
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FIG. 14: (Color online) Color coded local density of states 
pi{uo) as a function of frequency uj and lattice site index i 
for three different parameter sets: (a) U = 1.0 and sd = 
0.01E R , (b) U = 3.1 and s D = O.OIEr and (a) U = 1.0 and 
sd = 0.04Er. Parameters are v — 1.0 and the lattice size 
was 24 x 24. 



D. Finite temperature 

Here, we investigated a system of fermions with a semi- 
elliptic DOS at finite temperature. The spectral func- 
tions are plotted in Fig. [15] for two parameter sets. Panel 

(a) displays the arithmetically averaged spectral func- 
tions for speckle disorder strength sp = 0.05 and scat- 
tering length a = 150ao for various increasing temper- 
atures. At zero temperature, this parameter set would 
correspond to the Mott insulator. In Fig. [15] we note 
that with increasing temperature, the gap initially grows 
and the incoherent excitations reveal a significant redis- 
tribution of the spectral weight, which is shifted away 
from the Fermi level. This is not observed in the homo- 
geneous case, where the spectral transfer is weak. Panel 

(b) shows the corresponding behavior for stronger disor- 
der, namely sjj =0.1 and a = 150ao, which corresponds 
to the disordered strongly correlated metal at zero tem- 
perature. With increasing temperature, the coherent low 
energy excitations are reduced, and for /c#T = 0.05 the 
system enters a Mott insulating state. 

In Fig. [l6j the evolution of the arithmetically aver- 
aged spectral function with increasing temperature is 
displayed for two values of the speckle field strength 
(sd = O.IEr in panel (a) and sd = 0.2Er in panel (b)) 
and three different temperatures: /c#T = 0, ksT = 0.01 
and ksT = 0.05. We note the reduction of the spectral 
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FIG. 15: (Color online) Evolution of the arithmetically aver- 
aged spectral function with increasing temperature for two 
parameter sets: (a) sd = 0.05Er, a = 150ao and (b) 
sd — O.IEr, a = 150ao. Parameters are v — 1.0 and 
s L = 10E R . 



weight at the Fermi level due to finite temperatures. This 
leads to an enlargement of the Mott insulating phase. 
However, as can be seen in panel (b) of Fig. 16 , a metal 
without a resonant peak at the Fermi level is stabilized 
a higher disorder strength. In this respect, our central 
finding at T = 0, that the Mott and the Anderson-Mott 
insulators are not continuously connected in presence of 
the speckle disorder, also holds at finite temperature. 
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FIG. 16: (Color online) Evolution of the arithmetically 
averaged spectral function with increasing s-wave scatter- 
ing length (a = 50ao, 100ao, 150ao, 200ao) for two disorder 
strengths (a) sd = 0.1 Er and (b) sd — 0.2Er. For each 
parameter set the spectral function is compared at three 
different temperatures, namely kbT = (black solid line), 
kbT = 0.01 (red dashed line) and kbT = 0.05 (green dash- 
dotted line). Parameters are v — 1.0 and sl = IOEr. 



gap was found, which should be investigated in detail in 
the future. Finally, we investigated the high-dimensional 
system at finite temperature, where the Mott insulating 
region is enlarged, but the separation of the Mott and 
the Anderson-Mott insulators persists. 



SUMMARY AND OUTLOOK 
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In conclusion, we have investigated a gas of ultra- 
cold fermions in an optical lattices subjected to an ad- 
ditional speckle disorder field using statistical DMFT 
and RDMFT. The presented DMFT schemes include off- 
diagonal hopping disorder and allow for a systematic in- 
clusion of correlations between the difference in neighbor- 
ing on-site energies and the hopping amplitude, as well 
as between the on-site energy and the local interaction 
strength, which are non-trivial in realistic experiments 
and have been calculated by Zhou and CeperleyP^ 

The complete paramagnetic ground state phase dia- 
gram has been determined. It consists of a disordered 
metallic phase, as well as Mott insulating and Anderson- 
Mott insulating phases. A strong suppression of the 
correlation-induced metal insulator-transition by disor- 
der is observed and a finite metallic phase is found, even 
in the strongly interacting regime. Hence, the Mott 
and Anderson-Mott insulators are not continuously con- 
nected, in contrast to the prediction by typical medium 
theory in the presence of bounded box disorder. 

We have also investigated speckle-disordered fermions 
on a square lattice by means of RDMFT. We find that 
our main results for the case of high-dimensional systems 
also hold in two spatial dimensions. Moreover, a pseudo- 
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Appendix A: Sampling the Hubbard parameter 
distributions 



The various Hubbard parameters in an optical lattice 
exposed to speckle disorder with a realistic autocorrela- 
tion length of 1/im, such as realized in the experiments 
in the DeMarco group for instanc e,^ w ere calculated by 
S. Q. Zhou and D. M. Ceperle}ff^ at fixed disorder 
strength sd = ^Er and lattice intensity sl = 142£#. 
Imaginary time evolution was used, calculating the pa- 
rameters J7, e and t for 1222 disorder realizations on a 
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three-dimensional 6x6x6 lattice. Based on these data, 
PDFs were constructed, which are sampled and scaled 
appropriately in our calculations for the various param- 
eter regimes of interest. 

This scaling behavior will shortly be elucidated in the 
following: The most likely value and lower bound for 
the on-site energy corresponds to regions of low intensity, 
where we set e = 0, thus recovering the usual energy ref- 
erence point in the absence of the disorder speckle laser. 
It is found in Ref. ([T6]) that the disorder strength A (i.e. 
the standard deviation of p e (t)) scales linearly with the 
speckle intensity s^, as is to be expected in the perturba- 
tive limit. Thus, after an on-site energy is sampled from 
the distribution, it is simply scaled by multiplication with 
the respective sjj. An afflne shift is not required due to 
the vanishing lower bound discussed above. 

Since the on-site energy e and interaction parameter U 
are correlated, as can be seen from the histogram Fig. [3] 
(i.e. this two-dimensional PDF is not reconstructible as a 
product of two one-dimensional PDFs), a combined sam- 
pling from a two-dimensional distribution is used for e 
and U. This is equivalent to first sampling e from the 
marginal distribution p e {e) and subsequently U from the 
conditional PDF pu(U\e). Furthermore, the standard de- 
viation of the on-site interaction parameter U (i.e. of 
the marginal distribution pu(U) = J dep e jj(e,U)) scales 
linearly with s^, while the most probable value of U 
remains unaffected by a variation in sd up to a good 



approximation. Therefore, the variation in U in the two- 
dimensional distribution p e ,t/(e, U) is scaled by sp, while 
the most likely value is set to the value of U determined 
from a band structure calculation of the Wannier func- 
tions in the pure case without disorder, as performed in 
Ref. (44) for instance. 

In contrast to the local interaction parameter C7, the 
nearest neighbor tunneling amplitude tij exhibits only 
very weak correlation with the respective on-site energy, 
while the correlation with the energy difference between 
the two sites Ae = t{ — ej is significant. Hence, the dis- 
tribution for Ae and t cannot be sampled independently 
(i.e. £>Ae,t(Ae,£) ^ p^ e (Ae) pt(t)) and a conditional dis- 
tribution function p t (t\Ae) for t, given a fixed value of 
Ae, was constructed from the data irP4 For 200 discrete 
values of Ae a histogram approximating p t (t\Ae) was ex- 
tracted, approximating the PDF. This PDF is integrated 
with respect to t, yielding the conditional PDF F t (t\Ae) 
and subsequently normalized for each fixed Ae, such that 
lim^oo Ft(t\Ae) = 1. To randomly sample values in a 
numerically efficient manner from p t (t\Ae), the condi- 
tional cumulative PDF F t (t\Ae) is inverted with respect 
to t on a linearly interpolated grid on [0, 1], consisting of 
800 points. Given a fixed Ae, a randomly drawn value of 
the inverted cumulative distribution is thus distributed 
according to the initial conditional PDF p t {t\Ae), leading 
to the sought-after sampling routine. 
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